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Abstract 

Holas and March (Phys. Rev. A51, 2040 (1995)) wrote the gradient of the one-body 

potential V(r) in terms of low-order derivatives of the idempotent Dirac density matrix 
built from a single Slater determinant of Kohn-Sham orbitals. Here, this is first com- 
bined with the study of Dawson and March (J. Chem. Phys. 81, 5850 (1984)) to express 
the single-particle kinetic energy density of the Be atom ground-state in terms of both 
the electron density n(r) and potential V{r). While this is the more compact formula- 
tion, we then, by removing V{r), demonstrate that the ratio t{r)/n{r) depends, though 
non-locally, only on the single variable n'{r)/n{r), no high-order gradients entering for 
the spherical Be atom. 
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In recent work, we have been focussing on approximations to the exchange energy density in 
DFT [Ij. In this area, the long term aim is to express the single-particle kinetic energy T5 
and the exchange energy (plus later correlation energy) as functional of the ground-state 
electron density f2] . A recent review on such an orbital- free DFT has been presented by one 
of us in this Journal. [3] 

Here, we focus primarily on the kinetic energy Ts which we shall write in the form 



where t{r) is evidently the single-particle kinetic energy density. Throughout this Letter, 
we shall use the positive-definite form of t{r), corresponding to (V^E^)^ in terms of the one- 
electron wavefunctions. 

The general problem of finding Ts[n\ remains formidable. Therefore our focus here is on 
the example of the spherically symmetric Be atom ground state, with the configuration 
(ls)^(2s)^. Then, we know from the early study of Dawson and March [4j that the kinetic 
energy density t{r) entering in eqn. ([1]) can be usefully separated into two parts : 



where tiy(r) is the von Weizsacker ^ kinetic energy density given by 



So from eqns. (|2]) and ([3]) it is already clear that the sought-after functional form of Ts[n\ 
must involve both the density n and the gradient-dependent ratio \/n/n. 
Returning to Dawson and March wrote the Is and 2s Slater-Kohn-Sham [6l [7] orbitals 
in the forms 




(1) 



t{r) = tw{r) + - n{r) e'^{r) 



(2) 



22\E'is(r) = n^cosOir) 



(4) 
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and 

2^\l/2s(r) = n^sin6{r) (5) 

for normalised \l/'s. But, returning to eqn. (|2]), the most interesting result is a non-linear 
pendulum-like differential equation for the phase 9{r): namely 

Vn 

V^9{r) + .Ve{r) + Asm2^(r) = (6) 

It is immediately clear from eqn. ([6]) that the only input required to determine 9{r) in the Be 
atom is the gradient ratio Vn/n already mentioned in connection with the von Weizsacker 
component twif) in eqn. (|3]). Then the eigenvalue A, having dimensions of energy, and 
the phase 9{r) must be determined by integrating (usually numerically!) the non-linear 
(ordinary, for spherical atoms like Be or 4-electron atomic ions) differential eqn. dHJ. 
Such integration presumably will mean that while eqn. ([6]) shows already that 6* is a functional 
of Vn/n, i.e. 

e = 9[Vn/n] (7) 

the dependence of on r will not be determined by the 'local' quantity Vn(r)/n(r) at the 
point r at which 6 is to be evaluated. Thus, eqn. ([7]) is to be interpreted as a 'non-local' 
functional of the variable Vn/n. 

We turn at this point to introduce the one-body potential V(r). This enters, most basically 
in DFT [2j, through the chemical potential (//) equation 

Of course, in current applications of DFT, the functional derivative 6Ts/Sn{r) is bypassed 
by utilizing the Slater-Kohn-Sham (SKS) orbitals generated by the potential V^(r)^. Un- 
fortunately, this potential has an additive component from exchange and correlation, Vxcir), 
which is currently not known. 



For Be, we shall also bypass eqn. because of the difficulty associated with forming 
6Ts/6n{r) after inserting the formally exact eqn. (Q. What we appeal to instead is the 
density matrix formulation of DFT given by Holas and MarchjS] and used very recently by 
March and Nagy[n] in conjunction with the concept of the Pauli potential [lU]. From eqn. ([T]) 
of [8] we can write in spherical symmetry the force equation corresponding to the one-body 
potential V(r) : 

dVjr) _ l (|:)[V^n(r)] ^ z.jr) 
dr 4 n{r) n{r) 

The quantity Zs{r) entering eqn. is the single-particle (s) limit of the many-electron vector 

field z{r) defined in [S] and written explicitly in eqns. (2) and (3) of p. 

Utilizing the idea of the Pauli potential, in [9] it is shown, but now specifically for the Be 

atom, that 

. 1^ 2^ (10) 
r or 

where t{r) is given in eqn. ([2]) above. Thus, at the expense of working in eqn. ([9]) with the 
force, rather than V{r) in the chemical potential result ([8]), we have bypassed the need for 
functional differentiation of Tg [n] : of course the result ( 1TU|) being specific for the two-level 
Be atom. [TT] 

We must note here that in the very early work on the virial theorem in one dimension, [TT] 
it was shown that a differential virial form existed, namely 

dt 1 , .dV(x) 1 s 

It was clear already from the one-dimensional form ffTTj) that t{x) was very compactly ex- 
pressed in terms of both n{x) and its derivatives, plus the force —dV{x)/dx. Eqn. (fTUj) . 
when inserted into eqn. ([9]), leads to a related conclusion for the Be atom. 
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However, to achieve the prime object of this Brief Report, we want now, despite some 
mathematical comphcations, to avoid the force term —dV{r)/dr on the LHS of eqn. ([9]). To 
do so, let us return to eqn. |2]and define a kinetic energy difference A(r) by 

^"^'^'^ -'^^'^"^ = A(r) = e'\r) (12) 
n(r) 

Taking the square root of eqn. f|T2l) and noting that 6' (r) for Be is negative while A(r) is 
greater than zero, we obtain by integration that 

A^s)ds (13) 

This is the point at which to return to the non-linear pendulum-like eqn. (E]). Written for 
the spherically symmetric Be atom, we then find 

9" (r) + {q + -}9' (r) + A sin 29(r) = (14) 
r 

where q is the gradient quotient n {r)/n{r). 
Multiplying eqn. ( IT^ throughout by 9' then yields 



9\r)9"{r) + {q+ -}9"'r + X9\r)siYi29{r)9\r) = (15) 



2. . 



r ' 



The first two terms in eqn. (1151) can be written in terms of A(r) defined in eqn. (1121) . and 
its first derivative A'r, while the third term involves 9'r = — A^(r) as well as 9{r) given in 
eqn. f[T^ . The resulting equation then reads 

a' ( 2 1 /"^ 1 

+ {g + -}A'r + AA^ sin(-2 / A^ {s)ds) (16) 

Evidently, with input of g = n/n and also A, eqn. f lTB]) can be formally integrated to yield 
A(r). From eqn. (fT2|) this will yield the functional form 

t{r) 



n{r) 



F[g(r);A] (17) 



since from eqn. in spherical symmetry we have 

q\r). (18) 



tw{r) 2, 



n{r) 

Again, as in the formal functional form in eqn. ([7j), makes it clear that t{r)/n{r) depends 
only on the r space variable g(r) = n'(r)/n(r), which of course is a major simplification. 
It is worthy of note that, by a further appropriate differentiation using eqn. (|T6|) . A can be 
eliminated, but the mathematical detail proliferates and we shall not record it therefore. 
While this Letter has focussed on the kinetic energy density t{r), we conclude by noting that 
the exchange-only energy density ex{r) of DFT is customarily defined, following Dirac,[?] in 
terms of an idempotent first-order density matrix 'y{f,r). For Be, this has the form in terms 
of \E'is('") and\l'2s(r) written in eqns. (jl]) and (jS]) as 

7(r, r") = 2[^i,(r)^i,(r') + ^2s(r)^2.(r']) (19) 

But this is well known when using the forms (jlj) and in eqn. fll9p as 

7(7% r) = n{r)^n{r 

Evidently since the Dirac[T2] result for ex{r) in terms of 7(r, f^) is 

..(r) = 4/I^A^£*-' (21) 

4 J \r — r \ 

all we have done above on 6{r) for the Be atom is relevant to the calculation of the exchange 
energy density £xif) in terms of n(r) and the quotient g(r) = n {r)/n{r). 



' ' '''^cos[e{r)-9{r)] (20) 
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